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Abstract: 

A new Monte Carlo method for computing thermodynamical properties of very large polyelectrolytes 
is presented. It is based on a renormalization group relating the original polymer to a smaller system, 
where in addition to the naively rescaled forces, a corrective nearest-neighbor interaction originating 
from the short distance Coulomb cutoff is introduced. The method is derived for low T but is in the 
unscreened case valid for all T. Large polymers with N monomers are emulated by Monte Carlo 
calculations on smaller systems, K = N/Q. The computational gain of the method is Q^. It is 
explored with emphasis on room temperature. Results for 7\r=10000 are presented. 
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The thermodynamics of polyelectrolytes consisting of linear chains of monomers, with covalent 
harmonic bonding forces and Coulomb interactions, have been extensively studied with Monte Carlo 
Methods. Recently, high statistics results have emerged for relatively long (7\r=2048) chains [1]. A 
reliable measurement of e.g. the end-to-end distance r^e is very CPU time consuming, even when 
employing state-of-the-art updating methods [2, 3]. 

In order to circumvent the computational bottleneck, and to facilitate estimates of free energies, a 
variational scheme was proposed [4, 1], relying upon an Ansatz in the form of a generic Gaussian 
distribution. For the unscreened case the variational method overestimates r^e by 5-6 % at room 
temperature whereas for the screened case the approximation is less accurate. 

In this letter we develop an alternative Monte Carlo (MC) approach, where the chains are "coarse 
grained" in the spirit of the renormalization group. This is accomplished by using an energy ex- 
pression for blocked monomers with appropriately renormalized harmonic and Coulomb terms and 
an additional corrective nearest-neighbor interaction. Superficially, this is reminiscent of the well- 
known blob-picture of de Gennes [6], but the idea is quite different. 



The approach is designed for low T. Indeed, at least in the absence of screening, turns out to be 
low for large N [5]; thus, the approach should give good results precisely where it is needed. Our 
blocking scheme is derived in two ways, using (i) a real space renormalization philosophy, and (ii) 
a cut-off continuum formulation. 



In terms of dimensionless quantities, the energy of a Coulomb chain is given by [1] 

E = Ea + Ec = \Y.vl^, + Y.Y.^^'-^^) (1) 

« = 1 i jyi 

V{r) = —- (2) 

with a Boltzmann distribution oc exp(— iZ/T). Here, Vij = jx^— Xj| with as the monomer positions, 
while K is the inverse Debye screening length. With a length unit (defined by the unscreened N=2 
equilibrium distance) of 6A, room temperature (290K) corresponds to = 0.837808, and O.OIM, 
O.IM and l.OM salt solution correspond to k = 0.1992, 0.630 and 1.992, respectively [1]. 



Naively, one would coarse-grain the chain into blocks of Q monomers, and represent each of them by 
a single effective position carrying its entire charge. The remaining position variables can then be 
integrated out (being Gaussian), leaving a system of K = N/Q effective monomers, each carrying a 
charge Q. They would be connected by harmonic bonds being a factor Q weaker than the original 
ones. 



This would be fine, were it not for the divergent short-distance behavior of the Coulomb interaction. 
This mainly affects neighboring blocks, where the naive effective interaction underestimates the 
repulsion. 



This shortcoming can be remedied by a corrective interaction term between neighboring effective 
monomers, and we are lead to consider the following blocked energy: 

= ^ E + 1{N-K)T + Q'J2J2 ^(^^^O + E nri,i+i) (3) 

i = l i jyi « = 1 
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where Vij are distances between blocked monomers. The different terms represent, in order of 
appearance: (1) the naive effective bond energy, (2) a correction to the bond energy, to account for 
the bond energy of the N-K eliminated position variables, (3) the naive effective interaction, and 
(4) a corrective nearest-neighbour interaction. 

The correction term is determined as follows. In the low temperature limit the chain is a straight line 
with a slowly varying nearest-neighbour distance, which we approximate by a constant a. For the 
blocked system, the corresponding nearest-neighbour distance is then h = Qa. In this approximation, 
the true Coulomb energy is given by 

N 

Un = - (4) 

1=1 

while the blocked interaction energy (including the correction term) reads 

K 



Uk = J2^K - l)V{lh) + KW{h) 



1=1 



The objective is to choose W{h) such that ?7jv and Uk become identical. For the Coulomb potential 
of eq. (1) one obtains for Ui^: 



and similarly for Uk '■ 

UK^-^log{l-e-'''') + KW{b) (7) 

where only the leading term in the sums have been included. Equating ?7jv with Uk for b = Qa 
yields 

^(^) = - V N " '""^'^) " " 

In the absence of screening (k=0) this becomes 

r 

representing an enhancement of the nearest-neighbour Coulomb interaction. 

The results above are also transparent in a continuum formulation, where the ultraviolet divergence 
of the Coulomb interaction is regularized via an TV-dependent cutoff, as implied by the discrete 
nature of the original chain. The VF-term then corrects for the shift in the cutoff. We illustrate this 
in the unscreened case only. 

A continuum formulation of the chain is obtained by replacing the discrete index i by a continuous 
parameter a = i/N, such that = x(cr). Then the energy of eq. (1) (in the unscreened case) 
transforms into: 

2^70 ^ ' Jo J.+ l/N |X(<7')-X(<7)| ^ ' 
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Were it not for the short-distance cutoff {1/N) in a, needed to regularize the Coulomb potential, all 
the 7\r-dependence could be collected in a global factor N , by making the rescaling x(cr) — > N'x[a). 
This corresponds to a mere rescaling of the temperature, and we would have the naive scaling 

xr""' <^ Q-i"'"' (11) 

The 7\r-dependent cutoff introduces a scale-breaking, though, and for K < N the last term of eq. 
(10) can be rewritten as 



N2 r r dada' 

Jo Ja+lIN |x(0-') - X(C7 

Jo ia+l/if |x(0-') - X(C7)| Jo Ja+1 



/N 

+1/^ dada' 



In the low-temperature regime, where the chain is more or less linear, we can use the approximation 

x((t') - x((t) ^ {a' - a)k{a) (13) 
and obtain for the last term of eq. (12) 



r+^'^ dada' ^ ^2 r^^'^ dada' 

Jo Ja+ilN |x(o-') - x(c7)| ^ Jo y^+i/jv {a'-a)\i{a)\ 

Jo m\Jl/N ^' ' 'Jo m\ 

which is precisely the continuum version of the last term of eq. (3), using W of eq. (9). 

The correction to the naive blocked energy function was derived for low temperatures, and the 
method, as will be demonstrated below, indeed gives excellent results at low T, both for ree and 
Ec, for unscreened as well as screened Coulomb interactions. 

At higher temperatures, the polymer will be less linear, and it is not a priori certain, that the 
method will work there. For the unscreened case in the high T limit, the exact results for ree and 
Ec for large N can be expanded in a series in 1/T, yielding [5] 

<rle> ^ 3(7V-1)T+-^-^ (15) 
4 f2N^I'^ 

< Ec > : 



3 V TT Tl/2 

while Eg can be determined from the virial identity [1] 2Eg = Ec + 3{N — 1)T. We note, that 
the expressions are consistent with the naive scaling behavior of eq. (11), in that both ree/N and 
Ec/N are functions oiT/N only for large N . The relevant small parameter of the high-T expansion 
is obviously N/T. For the blocked system of K effective monomers, the corresponding results are 
identical to the order shown. Thus, in the unscreened case the high T performance is under control. 

In the presence of screening, the situation is somewhat different. The correct and blocked high-T 
expansions for Ec differ, and it is actually the contribution from the correction term that dominates 
[5]. Hence, in the screened case the blocking method is reliable only at low T. 
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N 


Q=8 


Q=4 


Q=2 


Q=4/3 


40 


1.10 


1.04 


1.01 


1.00 


80 


1.06 


1.03 


1.01 


1.00 


160 


1.05 


1.02 


1.00 


1.00 



Table 1: ree{K)/ree{N) as a function of N and Q = N/K at T=Tr for unscreened chains. 

We have made extensive numerical evaluations of the blocking approach. All MC runs, both stan- 
dard and blocked, were performed using the pivot algorithm [2]. Most results are based on 10^ 
thermalization sweeps and 10^ measured configurations. Some results were taken from ref. [1]. 
Initially exploratory comparisons were made in order to find Q, T and k ranges where the blocked 
approach is trustworthy. It turns out that realistic values for T (T^) for unscreened chains are within 
the domain of application even with very high Q- values. As expected, in the presence of screening 
the blocking approach breaks down when the screening length becomes smaller than the blocked 
resolution. The results are nevertheless encouraging since even with strong screening (k=1.992), 
Q < 10 is feasible. 

When investigating the power of the method, the relevant parameter is Q = N/K. However, in 
addition one expects finite- 7\r corrections due to end-point effects. This is illustrated in table 1, 
where the ratios between ree computed with blocked and standard method, ree{K) and ree{N) 
respectively, are shown for T=Tr and k=0.0. As can be seen from table 1 finite N effects diminish 
as N increases and should not pose problems for large system sizes, 1000 — 10000 monomers, at 
which the method is aimed. 

We next compare the power of our blocking approach on ree and Ec for T = Tr and k=0.0 and 
0.630 respectively. The results for ree can be found in figs. 1 and 2. The corresponding energies, 
which very well approximate the "true" MC values, are found in table 2. 



N 


80 


160 


320 


512 


1024 


2048 


K=0.0 Ec (K) 
(N) 


2.46 
2.46 


2.76 
2.76 


3.04 
3.04 


3.22 
3.22 


3.48 
3.48 


3.72 
3.73 


K=0.63 Ec (K) 
(N) 


0.47 
0.47 


0.50 
0.48 


0.54 
0.48 


0.56 
0.48 


0.58 
0.48 





Table 2: Average internal Coulombic energies per monomer for unscreened and screened Coulomb 
potentials. N and K stands for full MC and blocked with K = 50 respectively. 

As can be seen from fig. 1 and table 2 the ree and Ec are amazingly well described by the blocked 
approach for the unscreened case. The value of ree is within 10% from the true value even for 
Q 100. For even larger N the situation should be even better with diminishing finite N effects. 
Indeed, using the blocked method we have estimated ree and Ec for a 7\r=10000 unscreened chain 
with K=250 and find ree/N = 2.04 and Ec/N = 4.25 The relative errors should be less than those 
quoted for K = 50 and N = 2048 (3.6% for ree and 0.3% for Ec)- The computational speedup here 
is substantial, growing as Q^. 
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Figure 1: r^elN as a function of N for an unscreened Coulomb chain. Full and dashed lines denote 
K=2Q and K=bQ blocked MC runs including logarithmic correction respectively. Also shown as 
dotted lines are the naive blocked results for K=2Q. The squares denote results from original MC 
runs with N degrees of freedom. 
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Figure 2: r^elN as a function of N for a k=0.630 screened Coulomb chain. Same notation as in 
fig. 1. 



For the screened case (see fig. 2) the approximation is less accurate but still good enough to be of 
practical relevance. In table 3 detailed comparisons are made for (5=4 and 8 for different screening 
strengths. As can be seen there, the blocking gives a reasonable approximation even for Q=8, giving 
a speedup factor of 512. 
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K 


r-ee(N) 


Q 


r-ee(K) 


% diff 


Q 


r-ee(K) 


% diff 


0.1992 


185 


4 


186 


0.6 


8 


187 


1.1 


0.630 


104 


4 


105 


1.0 


8 


110 


5.8 


1.992 


61 


4 


62 


1.6 


8 


68 


11.4 



Table 3: Comparison of T"ee(K,N) and ree{N,N) for N=512 screened chains with (5=4 and 8. 

The approximation is of course not relevant for quantities more local than the blocking allows for. 

In summary, we have developed an efficient Monte Carlo blocking scheme that allows for estimating 
certain thermodynamical quantities of polyelectrolytes with very substantial computational speedup 
factors for temperatures around T^.. The method is excellent for unscreened chains, but even with 
screening non-neglible savings result. Results for chain lengths hitherto never probed are reported. 
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